Psi Function

public function Psi(psic, theta, thetas, thetar, psdi) result(h)

Compute matric potential (m)

Arguments

Type IntentOptional Attributes Name
real(kind=float), intent(in) :: psic

bubbling pressure (m)
real(kind=float), intent(in) :: theta

volumetric water content (m3/m3)
real(kind=float), intent(in) :: thetas

saturated volumetric water content (m3/m3)
real(kind=float), intent(in) :: thetar

residual volumetric water content (m3/m3)
real(kind=float), intent(in) :: psdi

Brooks & Corey pore size distribution index (-)

Return Value real(kind=float)


Source Code

FUNCTION Psi &
!
 (psic, theta, thetas, thetar, psdi) &
!
 RESULT (h)

IMPLICIT NONE

!Arguments with intent in
REAL (KIND = float), INTENT(IN) :: psic !!bubbling pressure (m)
REAL (KIND = float), INTENT(IN) :: theta !!volumetric water content (m3/m3)
REAL (KIND = float), INTENT(IN) :: thetas !!saturated volumetric water content (m3/m3)
REAL (KIND = float), INTENT(IN) :: thetar !!residual volumetric water content (m3/m3)
REAL (KIND = float), INTENT(IN) :: psdi !!Brooks & Corey pore size distribution index (-)


!local declarations:
REAL (KIND = float) :: h

!-------------------------end of declarations----------------------------------


 h = psic / ( (theta - thetar) / (thetas - thetar) ) ** (1./psdi)

RETURN

 END FUNCTION Psi